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Abstract 

We propose a model for protein folding in vivo based on a Brownian-ratchet mechanism in 
the multidimensional energy landscape space. The device is able to produce directed transport 
^ taking advantage of the assumed intrinsic asymmetric properties of the proteins and employing 

' (— l ' the consumption of energy provided by an external source. Through such a directed transport 

^ phenomenon, the polypeptide finds the native state starting from any initial state in the energy 

^ landscape with great efficacy and robustness, even in the presence of different type of obstacles. 

Q This model solves Levinthal's paradox without requiring biased transition probabilities but at the 

^ expense of opening the system to an external field. 

> 

00 
(N 

o 
a^ 
o 

> 
X 



1 



I. INTRODUCTION 



Following the classical work of Anfinsen [H [2] on protein folding, Levinthal jS] argued 
that if the protein multidimensional free energy landscape had a "golf course" like shape, 
with the hole being the native state, it would take essentially forever for a protein to fold 
correctly, provided the search were performed randomly. This is known as Levinthal's para- 
dox, since proteins in vivo fold very fast and with high efficacy. However, Zwanzig et al 
O |6] demonstrated that, within a random search, a small bias for configurations closer to 
the native state would be enough to solve the paradox. Alternatively, to explain the efficacy 
of folding, the concept of a funneled energy landscape was developed [3 El 13 1101 IHl [12] , 
with the further possibility of folding pathways. In such an scenario the process is driven 
mainly by free energy differences, that is, by "falling" through the funnel into the state of 
lowest free energy. This idea has been central in the study of protein folding. In this article, 
we propose a further alternative way out to the so-called Levinthal paradox, which does not 
necessarily needs a funnel structure nor introduces ad-hoc biased probabilities. We appeal 
to a Brownian motor or ratchet-like mechanism for the process of protein folding. This is 
a non-equilibrium process that requires, ffist, a free energy landscape with a ratchet-like 
asymmetry, and second, the consumption of energy from an external source [13]. We shall 
discuss below how these two essential ingredients may be justified. 

The ratchet mechanism to produce directional current or motion has become a paradigm 
of the interplay of nonlinear phenomena with Brownian motion and unbiased external noise 
in mesoscopic systems [TH UHl UHl |T7] . We shall use the term "ratchet" for any conserva- 
tive mechanical system that has an ingrained spatial asymmetry; a simple one-dimensional 
example is a particle in a saw-tooth like potential. Due to the potentiality of these systems 
to describe mesoscopic systems such as biological ones [151 IIHl [HI [H] , one should consider 
the ratchet to be immersed in a thermal bath. This adds dissipative and stochastic thermal 
forces that are related to each other through the fluctuation-dissipation theorem [T^l I2U] . 
As it was clearly pointed out by Feynman in his Lectures |21], the Second Law of Ther- 
modynamics implies that a ratchet under the above conditions cannot generate current nor 
directional motion. Nevertheless, it is by now very well established that in order to obtain 
directional motion, the presence of an external unbiased time dependent force is necessary. 
This external force can be deterministic, a so-called "rocking" ratchet [E], or stochastic in 
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nature [22]. In any case, any of those forces should be of zero average in time such that the 

current obtained is not a trivial consequence of the bias of those forces, but of their inter- 
play with the ratchet potential. In the present article, since the external interaction of the 
protein during the folding process is supposed to yield energy through the consumption of 
high energy phosphate bonds, as adenosine-triphosphate (ATP) or guanosine-triphosphate 
(GTP), we assume an external, unbiased, stochastic force. 

We shall show that, understanding by protein folding the search and finding of a target 
point or small section in the energy landscape, the ratchet mechanism is extremely effective 
and robust. As mentioned above, an interesting consequence of this process is that the 
landscape does not require a funnel structure [71 ITO] . 

II. THE MODEL 

The model consists of an overdamped Brownian "particle" moving in a multidimensional 
configuration space, the coordinates representing the effective degrees of freedom of a protein 
[231 . The environment enters as a dissipative thermal bath and the full folding process is 
driven by an unbiased external source. 

Our central proposal is that the regularity and intrinsic asymmetries found in polypep- 
tides, such as a helices and /5 sheets, may arise from or give rise to an asymmetry at the 
level of the energy landscape. This feature is introduced through an asymmetric ratchet-like, 
conservative potential. The potential may be flat on average but it is provided with a small, 
deep and wide enough well, whose local minimum vicinity represents the native state of the 
protein. On top of the average potential there is a ratchet-like asymmetry that makes the 
potential towards the native state different than in the direction moving away from it. Figure 
[l] is a two-dimensional "free energy" landscape representation of the potential. Although 
the precise form of the landscape is essential for a given protein, here we shall use a very 
simple form for the configuration potential in order to show the virtues of the model. In d 
dimensions the ratchet-like potential with a native state is given by. 



where x = Xi,X2, ■■■,Xd are the configuration coordinates; r = {xf + xl + ... + x'^)^; and a, 
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FIG. 1: (Color online) Two-dimensional representation of the protein "free energy" landscape, 
Eq.Q, with A = 1, Vb = 1, Fi = 50, a = 0.4, b = 0.1 and c = 10. This landscape is flat on 
the average and has a ratchet-like asymmetry away and towards the native state, represented by a 
deep and wide enough minimum of the potential. 

b, c, Vq, V\ and A are positive real numbers that determine the structure of the potential. 
In the figure captions we specify their particular values. The native state is arbitrarily set 
at r = and it is a well given by the last term in Vix). The first three terms give rise to a 
periodic ratchet potential with a radial asymmetry. 

The above landscape potential can be make extremely more complicated and still work 
in the same way as the stylized model. That is, nonperiodic ratchets with wells, barriers 
or other structures can be introduced in order to approximate the potential to the real 
conditions of proteins, without affecting the essence of the model. This will be shown below. 

Proteins are typically embedded in a viscous environment at a fixed temperature. This 
is represented in the model by a thermal bath exerting both a dissipative and a white, 
Gaussian stochastic force. We believe this does not require further justification. 

With the last two ingredients only, namely, with the ratchet potential and the thermal 
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bath, the "particle" will obey Levinthal paradox and would not find its way to the native 
state, particularly in a multidimensional space. The directed transport mechanism requires 
the presence of a time-dependent external source of energy acting on the particle in order to 
produce work, in the form of a directional current [2^. To represent a natural, non-designed 
process, this force must be unbiased. It turns out that the interplay of the asymmetry of 
the conservative ratchet potential with the external source of energy, yields the sough for 
directed motion. We shall demonstrate below how the introduction of this force makes the 
particle finds its way toward the native state and, on the average, in a very short time 
[13]. The efficacy of the process depends, of course, on the values of the different model 
parameters. 

Regarding the origin of the external energy source, we generally argue that processes 
in living organisms occur via non-equilibrium states with all sorts of gradients of different 
physical quantities, and that the maintenance of those gradients may be traced back to the 
consumption of energy. In a more specific way, during the protein synthesis, the action 
of the ribosome it is an important infiuence to the folding process, in the sense that this 
process is almost concluded in an average size protein by the time the ribosome releases its 
C-terminal end [26]. Furthermore, ribosomes are recognized complex molecular motors that 
consume energy in the form of GTP throughout the protein synthesis [27t l28t 129] . Although 
it seems it does not occur in all folding cases, molecular chaperones may represent another 
good example of an external "force" in the folding process since they consume ATP while 
assisting in the process [301 Ell 1321 1331 [31] • In a broader sense, we may also argue that of all 
the many processes occurring in the crowded environment of the cell, some of them consume 
energy and can be used by the protein during the folding process as the required external 
force to reach the native state. 

The mathematical model is, therefore, the overdamped Langevin equation for a Brownian 
particle embedded in a thermal bath moving in a multidimensional coordinate system, under 
the infiuence of a conservative potential V{x) and in the presence of a time-dependent 
external force Fi(t), 



with i = 1,2, . . . ,d. 7 is the friction coefficient of the bath; fi{t) is an stochastic Gaussian 
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thermal force, with zero mean and with a white spectrum, namely [T9]: 



{mfjit')) = 2^kT5it-t%„ (3) 

with i,j = 1,2, ... ,d. In this expression T is the temperature and k Boltzmann constant. 
As mentioned, {fi(t)) = for all i. Thus, f(t) represents unbiased white noise. The fact 
that the thermal force correlation is proportional to 7/cT is the statement of detailed balance 
between the thermal bath and the particle [191 ED] • 

As mentioned, the forces Fi{t) represent the action of agents external to the protein itself 
and the thermal bath, and its driving is caused by ATP consumption. We really do not know 
the precise form of this force. However, as far as producing directed transport, its precise 
form is not essential, as long as its correlation time is different from that of the bath . 
Therefore, we have decided to employ a stochastic force generated by an Ornstein-Uhlenbeck 
process [T9|, [20] given by, 

^^-f>.C.W, (4) 

for i = 1, . . . ,d and where the functions Q{t) represent independent stochastic variables with 
zero mean {Qit)) = 0, and Gaussian white noise, 

{mCAt')) = ^f,5it-t')S,,. (5) 

e 

The parameters Tg and fo determine the process and represent the correlation time and the 
fluctuations magnitude of the external force, respectively. As we pointed out above, it has 
been proved that Brownian motors also work with different kind of zero bias forces, such as 
sinusoidal or dichotomic ones ITS). 



III. RESULTS 

We solve Eqs.(|2|-([5]) numerically. We first fix the potential parameters, Vq, Vi, a, b, c and 
A, and then we vary 7, T, fo and Te such that the directed transport process is optimized. 
Due to numerical easiness, we study the evolution of a protein with six degrees of freedom. 
However, to show the robustness of the model we also present results for a protein with 100 
degrees of freedom. This may represent a real protein with several tens of aminoacids. 
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FIG. 2: Typical folding realization of a protein with 6 degrees of freedom. Each line refers to the 
evolution of a single coordinate Xi, with i = 1, . . . ,6. (a) The stochastic external force is included, 
(b) No external force is considered and no folding occurs. In this case, 7 = 1, kT = 0.1, /o = 3.0 
and Te = 0.5. The initial condition was 200 units in each coordinate for both cases. 

A. Folding in 6 dimensions 

Figures [2^ and [2]d show typical examples of the particle evolution in a 6 dimensional 
space, by following the value of each of the coordinates Xi for i = 1, . . . , 6. The behavior 
when the stochastic external force is present is shown in Fig. |2^, while Fig. [2]d corresponds 
to the case when no external force is considered (and the temperature is increased in order 
to compensate in magnitude the lack of the external force). An evident directed transport 
towards the native state is observed in the former case, in contrast with the random-walker 
behavior when no external force is added. The directed transport occurs (Fig. |2^) until 
the particle reaches the vicinity of the native state well. Then, a fluctuation takes place 
and the particle finally falls into the well. If the well is not wide enough the particle may 
stay for a long time around the minimum of the well. Once the particle is within the well, 
its position fluctuates depending on, the width of the well at its bottom, the magnitude of 
the two stochastic forces, and the dimensionality of the space. Typically, for spaces with 
higher dimensionality wider wells need to be considered in order to make the particle reach 
the native state. In Fig. |3| we show the evolution of the dimensionless radius position 
r' = (xi + X2 + ... + x^)5/ro, with Tq the initial condition. This figure corresponds to typical 
cases such as those analyzed above as well as the average for 100 realizations. For the 
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(d) r' without external field 
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FIG. 3: Time evolution of the dimensionless distance to the native state r' = (xf + 2 /ro, 

with To the initial condition, for a protein with 6 degrees of freedom. Single realizations and averages 
over 100 realizations are shown for both cases, with and without external forces. Same parameters 
as in Figj2| 

relevant case when the external force is included, cases (a) and (b) in the figure, we have 
found 100% efficacy of "folding" within the arbitrary time limit of 10000 units of time for the 
parameters considered. This time is very short even in CPU time scales. When the external 
force is absent, the particle is a true diffusing random walker and, insofar our numerical 
calculations, we have never seen even a single occurrence of finding the native state. This 
point serves to emphasize the role of the external force as the responsible one for delivering 
the energy to be further converted by the ratchet into useful work. The result of the work 
done is the folding of the protein. 

B. Robustness 

In order to demonstrate the robustness of the model two representative cases were stud- 
ied. The first consisted on increasing the degrees of freedom up to 100 to verify that the 
dimensionality of the space has no fundamental influence on the behavior of the Brownian 
motor. The second case considers the inclusion of obstacles into the energy landscape, such 
as potential barriers and wells, in order to evaluate the capacity of the folding process to 
overcome such expected features in more realistic scenarios. 
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FIG. 4: Typical folding realization of a protein with 100 degrees of freedom, with the external force 
present. Each color refers to the evolution of a single coordinate Xi, with i = 1, . . . , 100. Here, we 
used the same potential parameters as in Fig{T] except Vi = 10000 and c = 6000, while the other 
parameters are: 7 = 1, kT = 0.1, /o = 3.0 and Tg = 0.5. The initial condition is 300 units for each 
coordinate. 
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FIG. 5: Time evolution of the dimensionless distance to the native state r', see Figjs] averaged 
over 20 realizations for a protein with 100 degrees of freedom. Same parameters as in Figj4} 

1. 100 dimensions 

Figure |4] shows the evolution of the 100 state coordinates of a protein with 100 degrees of 
freedom. As before, the particle is under the influence of a potential such as that of Eq. ([T]) 
and the external field is included. Figure [5] shows the time dependence of the dimensionless 
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radius (r'), averaged over 20 realizations with the same conditions as the case shown in 
Fig. |4j In this case, the initial condition is 300 spatial units for each of the coordinates. 
We find that the overall behavior for 100 coordinates is essentially the same as for the 6- 
dimension case, showing an evident directed transport towards the native state with 100% 
of success reaching it. There are, as expected, quantitative differences as the dimensionality 
is increased. One is the increment in time to reach the native state, as seen in the figures. 
Other is the increase in size of the spatial fluctuations around the well. Both effects are a 
consequence of having much more space to explore in the directions orthogonal to the radial 
one, which is the unique direction where the ratchet asymmetry is present. To be precise, 
within the parameters used, in order for the particle to get into the native state well, it was 
necessary to increase the width of the well from 10 to 250 space units as increasing from 6 
to 100 dimensions. Further, in order to achieve small fluctuation amplitudes in the native 
state, see the zoom-in graph in Figjl], the width at the bottom of the well was kept narrow 
but its depth was increased. It should be also clear that the time to reach the native state 
depends on how far from it the particle is placed initially. At the model level this is a matter 
of CPU time only. 



2. 6 dimensions with obstacles 



The bidimensional representation of the energy landscape of a protein with two additional 
obstacles is shown in Fig. |6} The obstacles consist on, one, an additional Gaussian well (I) 
deeper and wider than the native state well (N), (implying in principle a more stable state) 
and, two, a Gaussian barrier between the two wells and surrounding the native state. The 
equation corresponding to this potential with obstacles in d dimensions is given by 



V{x) = Vo 



, 2TTr\ ( 47rr\ ( Gvrr 

sin — — + a sin — — + o sin 



\ ) \ \ ) \ \ 

Vie-'"'' - VbC-^''-''^^"'"^ + V^e-'"-''-, (6) 



with = {xi — x'^Y + {x2 — + ■ ■ ■ + {xd — x'^Y . The new parameters H, V^, and are 
the positive real numbers that determine the shape of the barrier and the well introduced, 
while Tf, and xlj = {x^, X2 , ■ ■ ■ , x'^) determine their positions, respectively. 

The behavior of a particle in a 6 dimensional space, under the influence of the potential 
with obstacles, Eq. ([6]), and subjected to the stochastic external field, is shown in Fig. [?} 
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FIG. 6: Bidimensional representation of the potential with obstacles, Eq. |6j Here, A = 1, Vq = 1, 
Vi = 50, a = 0.4, b = 0.1, c = 10, 14 = 18, K; = 60, Cb = 30, = 30, n = 250, and xf = 200 for 
i = 1 , . . . , (i. This conservative potential corresponds to the energy landscape of a protein with the 
native state (N) centered in the origin and two additional obstacles, a well (I) far away from the 
native state and a barrier between the two wells. 

It is important to mention that the initial condition was set at the minimum of the added 
well (I). Fig. [T] shows that the particle is able to escape from the first, deeper but wider 
well, then overcome the potential barrier and finally find the native state. Frequently, the 
obstacles are not overcome immediately, as shown in the radius position evolution of typical 
realizations denoted by (a), (b), (c) and (d). We see regions of constant distances revealing 
that the particle takes some time leaving the first well and climbing the barrier. However, 
the particle always reaches the native state within the arbitrary time limit of 15000 time 
units, in an average time longer than the case with no obstacles , as the average radius of 
100 realizations (e) illustrates. 

Therefore, the protein "folds" even if severe obstacles, such as barriers and potential 
wells, even deeper than that of the native state, are included in the energy landscape. This 
climbing capacity of the particle is another consequence of the delivery of work by the 
ratchet. The climbing capability does have a stall limit, thus, that depends on the local 
slope of the obstacle but not on its height or depth. 
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FIG. 7: Time evolution of the dimensionless distance to the native state r' , see Figjsj for a protein 
with 6 degrees of freedom in the potential with obstacles shown in FigjGj Cases (a)-(d) show typical 
single realizations and (e) is the average over 100 realizations. The initial state is 200 units in each 
coordinate and coincides with the minimum of the additional well (I). 

IV. CONCLUSIONS 

In this article we have revisited the ratchet or Brownian-motor mechanism extended to 
several dimensions in order to propose a model for protein folding. The ratchet mechanism 
has become a potential candidate for many biological processes, since, on the one hand it 
resembles a thermal engine at a mesoscopic level, but on the other hand, it does not appear 
to need a specific design such as the man-made engines. That is, it needs only of two essential 
ingredients: an intrinsic asymmetry and the input of external unbiased energy. The result 
is delivery of work in the form of directional motion. 

With this in mind several authors have used ratchets models in more obvious biological 
situations in which a transport process, of some kind, is present [TSl [ISl [13 HE]- Protein 
folding appears as a process with a directional motion towards the native state in the mul- 
tidimensional space of the variables that define the energy landscape. These facts persuade 
us to propose a model for the mechanism of protein folding based on a Brownian motor. 

Traditionally, this process has been thought to be driven by free energy differences leading 
towards a minimum, similarly to a chemical reaction [8], and thus the idea of a funnel in 
the free energy landscape of the protein is very appealing. The contribution of external 
agents are welcome in the theory, however, they are not required. Here, by following the 



12 



idea that Life processes consume energy to yield their products, we have speculated that 
protein folding is driven by an "engine" in which the protein itself is part of it. Therefore, 
the presence of external sources of energy, provided by the ribosome, molecular choperones 
or other elements which consume energy during the process, becomes a necessity for the 
model presented here. 

We have focused our work to the in vivo protein folding phenomenon since there are 
more evident possible sources for the required external energy. However, we think that the 
model presented can also be considered for the in vitro situation, where part of the bias and 
external energy necessary to obtain directional work to refold the proteins, may be provided 
by the laboratory procedures, such as changes of temperature, pH or concentration of urea 
[35| [36] , among others. 

On the other hand, at this moment the most difficult part to justify in our model is 
the asymmetric structure of the energy landscape. The best that we can say now is that 
the intrinsic asymmetric structure of proteins, which also presents certain periodicity, en- 
courages us to think that ratchet-like structures can be hidden into the energy landscape. 
Nevertheless, we do not have evidence neither pro nor con that a ratchet structure is present 
since this would have to be seen in the highly multidimensional energy landscape. 

By means of numerical solutions of the model, it was found that the present model in 6 
dimensions is able to produce directed transport in the mentioned space, through which the 
polypeptide is able to found the native state starting from an arbitrary initial state in the 
energy landscape and with a 100% success. This happens when the external source of energy 
is present, otherwise the finding of the native state is essentially never achieved. When the 
number of degrees of freedom is extended to 100, again, the process is performed with 100% 
efficacy, but due to the larger number of dimensions, the time needed to reach the native 
state increases approximately 6-fold with the parameters here used. On the other hand, 
the model demonstrates to be quite robust by showing that the protein is able to overcome 
severe obstacles in the energy landscape in its way to the target state. 

To summarize and to emphasize it again, the interesting and powerful side of this model 
is that it achieves its goal by taking advantage of two simple requirements only: the intrinsic 
asymmetric properties of the system and the external energy source. The numerical solutions 
showed in this work prove that Levinthal's paradox can be solved with a model which does 
not require ad hoc biased transition probabilities. In contrast, it requires the necessary 
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driving of an external energy source, which in our opinion is physically appealing since its 
presence is essential for Life to occur. 
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